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Abstract — The population synthesis method is used to study the possibility of explaining the appreciable 
fraction of the intergalactic type-la supernovae (SN la), 20t\l%, observed in galaxy clusters (Gal-Yam 
et at. 2003) when close white dwarf binaries merge in the cores of globular clusters. In a typical globular 
cluster, the number of merging double white dwarfs does not exceed ~ 10 -13 per year per average cluster 
star in the entire evolution time of the cluster, which is a factor of ~ 3 higher than that in a Milky- Way- type 
spiral galaxy. From 5 to 30% of the merging white dwarfs are dynamically expelled from the cluster with 
barycenter velocities up to 150 km s -1 . SN la explosions during the mergers of double white dwarfs in 
dense star clusters may account for ~ 1% of the total rate of thermonuclear supernovae in the central parts 
of galaxy clusters if the baryon mass fraction in such star clusters is ~ 0.3%. 



INTRODUCTION 

SN la (thermonuclear supernovae) are among the 
most important objects in modern astrophysics, be- 
cause the maximum of their optical light curves has a 
small spread. Therefore, they can be used 1 as a stan- 
dard candle in modern cosmology, and fundamental 
conclusions about the expansion kinematics of the 
Universe can be reached (Riess et al. 1998; Perl- 
mutter et al. 1999). The physical cause of an SN la 
explosion is the thermonuclear explosion of a white 
dwarf (WD) with a mass close to the Chandrasekhar 
limit (Hoyle and Fauler 1960). This idea is fully con- 
firmed by detailed numerical calculations of the evo- 
lution of a white dwarf in a binary system (see, e.g., 
Dunina-Barkovskata et al. 2001 ). The WD mass can 
increase during mass transfer between the compo- 
nents in a close binary, which determines the astro- 
physical sources of such supernovae — double WDs 
that merge through gravitational radiation (Iben and 
Tutukov 1984; Webbink 1984) or accreting WDs in 
semidetached or symbiotic binaries with a second 
nondegenerate component (Whelan and Iben 1973; 
Nomoto 1982). SN la are observed in galaxies of all 
morphological types (van den Berg et al. 2003), and 
their rate (in SNU units) depends weakly on the type 
of galaxy (see Capellaro et al. (1997) and more recent 
papers of Capellaro's group). 
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1 However, with certain reservations, up to 40% of SN la may 
be peculiar (Lee et al. 2001), plus a significant spread of 
theoretical light curves (Sorokina et al. 2000). 



At present, it is not completely clear which of these 
scenarios (or both) make the largest contribution 
to the SN la rate, and a considerable observational 
effort is being made to find the Galactic population 
of double WDs with a total mass close to or larger 
than the Chandrasekhar limit (e.g., the SPY project 
by Napiwotzki et al. 2003). We will adhere to the 
scenario of merging double WDs as the progenitors of 
type-la supernovae. The recent discovery of hydrogen 
emission lines in the spectrum of SN 2002ic (Hamyu 
et al. 2003) cannot rule out this formation channel 
of thermonuclear supernovae, because such SN la 
are atypical; their origin is the subject of debate (e.g., 
Livio and Riess 2003; Chugai and Yungelson 2003). 

A purposeful search for supernovae in Abell 
clusters of galaxies (the WOOTS project; Gal-Yam 
et al. 2003) has revealed that two of the seven 
SN la found (SN 1998fc in Abell 403 at z = 0.1 
and SN 2001al in Abell 2122/4 at z = 0.066)) have 
no obvious host galaxies and can be associated 
only with faint dwarf galaxies with Mr > —12™. 
In reality, in both cases, the supernovae are pro- 
jected onto the halos of the central cD galaxies in 
the corresponding clusters. However, the authors 
rejected the possibility that they belong to these 
galaxies because of the large radial-velocity difference 
(750-2000 km s -1 ) between the supernovae and the 
corresponding cD galaxies. 

Assuming a modified Schechter luminosity func- 
tion for the galaxies in the Virgo cluster (Trentham 
and Tully 2002), the luminosity fraction of the dwarf 
galaxies in the clusters under study should have been 
less than 0.3%. If the supernova rate is proportional 
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to luminosity, and if the fraction of the supernovae 
without obvious host galaxies is on the order of 10%, 
this leads to their ~ 30-fold excess among such dwarf 
galaxies compared to high-mass galaxies. Therefore, 
these supernovae are believed to have exploded out- 
side galaxies (intergalactic supernovae) as a result of 
the evolution of the intergalactic stellar population. It 
is concluded from these observations that the fraction 
of the intergalactic SN la is 201^%, in accordance 
with the theoretically expected fraction of the stars 
outside galaxies due to tidal interaction between clus- 
ter galaxies. The authors themselves (Gal-Yam et 
al. 2003) point out a low statistical significance of 
the result obtained, and such a high percentage of the 
extragalactic SN la should be confirmed by further 
observations. 

There is a potential channel for the formation of 
intergalactic supernovae through the dynamical ex- 
pulsion of close WD pairs from the system of galactic 
globular clusters (GCs). The stellar mass fraction in 
the observed GCs is ~ 0.1% of the visible mass of the 
galaxies. However, dynamical interactions between 
stars in dense cluster cores increase the merger rate 
of close pairs with white dwarfs. Some of them can 
be expelled from the clusters during tripple collisions, 
and the assumption that the rate of thermonuclear 
supernovae is proportional to luminosity for dense 
clusters becomes invalid. The numerical calculations 
of the WD evolution in open star clusters by Shara 
and Hurley (2002) confirmed the significant increase 
in the formation rate of close WD pairs through dy- 
namical interactions. 

The goal of our study is to quantitatively estimate 
the merger rate of double WDs in the dense cores 
of star clusters by using the method of population 
synthesis of the evolution of binary stars with a semi- 
analytical allowance made for their dynamical inter- 
actions with single stars in the cluster core. We show 
that the the merger rate of double WDs in GCs per 
average cluster star is a factor of ~ 3 higher than the 
merger rate of double WDs in spiral galaxies with 
continuous star formation. If ~ 0.3% of the baryon 
mass in the cluster centers is concentrated in virial- 
ized dense star clusters with masses of 10 5 — 10 8 M Q , 



then ~ 1% of SN la in the cluster centers can be 
naturally explained without additionally assuming the 
presence of 10% of the intergalactic stellar popula- 
tion. 



THE MODEL 

The Structure and Evolution of a Globular Cluster 

In our calculations, we used a Michie— King model 
(Michie 1963) to describe the structure of the globular 
cluster. The stellar population was broken down into 
subsystems (depending on their mass); their space 
density p a was fitted by a power law: 
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where r Ca and are, respectively, the radii of the 
core and the sphere within which half the mass of the 
a subsystem of stars is contained; and r t is the tidal 
radius of the cluster. Thus, p ha = p Ca { r h a /r Ca )~ 2 
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energy equipartition condition for the GC stars, we 
may assume that 
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is the rms space velocity, and p c is the central number 
density of the stars: 
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The stellar velocities can be described by the trun- 
cated Maxwellian distribution 




Ke^ j2 [e-P E -e~P E % E < E t , J < J C (E), 



E > E t , 



(3) 



where E and J are the energy and the angular mo- energy E. Stars with E > E t are assumed to have 
mentum (per unit mass), respectively; and J C (E) is 

the value of J for a star in a circular orbit with escaped from the cluster. 
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Detailed calculations of the evolution of the stan- 
dard model reveal a collapse of the inner core and an 
expansion of the outer layers. As follows from numer- 
ical calculations (see, e.g., Kim et al. 1998), the core 
contraction at late stages asymptotically approaches 
a power law: 

pS) ~ p c (0)(l - t/tcoii)" 1 - 2 , v m (t) (4) 
~v m (0)(l-t/t coll )-° A2 , 

where t is the current time, and t co]] is the collapse 
time of the cluster. 

As the GC core contracts, the central number 
density of the stars rapidly increases, which causes 
the formation rate and the binding energy of bi- 
nary stars to increase. Once the initial collapse has 
stopped, one might expect an expansion of the cluster 
core due to the tidal formation of binary stars. For 
some time, the newly forming binaries will be the 
energy source that maintains the expansion of the 
inner cluster regions. The expansion can be described 
by the following asymptotic formulas (see Kim et 
al. 1998): 

Pc (t)~t- 2 ,v m (t)~t-°- S2 . (5) 

We used several GC models that differed in central 
density p c and radius r c at the time of the collapse. In 
all models, the collapse time was taken to be t co \\ — 
7 x 10 9 yr, the central potential was Wo = -(3(f)(0) — 
9((3 = <p(r) is the gravitational potential), the 

total number of stars was 10 6 , and the fraction of the 
initially binary stars was assumed to be 30%. 

The last two parameters specify the normalization 
of our calculations. For the adopted Salpeter initial 
mass function with a minimum mass of the forming 
stars of O.1M , the mean stellar mass in the GC is 
O.38M , and the total mass of the GC of 10 6 stars 
is ~ 4 x 1O 5 M . Thus, the results presented below 
pertain to modeling the evolution of 3 x 10 5 binaries 
in a GC with a mass of 1.3 x 10 5 M Q . 

Here, it is pertinent to make a note on the chosen 
30% fraction of the initially binary stars in the GC. It 
follows from Hubble Space Telescope observations of 
GCs (Rubinstein and Bailyn 1997) that the fraction 
of the main-sequence binary stars in the GC with a 
collapsed core NGC 6752 may range from 15 to 38%. 
In NGC 288, this fraction is estimated to be in the 
range from 8 to 38% (Bellazzinni et al. 2002). A 
numerical analysis of the evolution of binaries in a GC 
(Fregeau et al. 2003) indicates that the fraction of the 
binaries with respect to their initial number after the 
core collapse ranges (depending on the cluster model) 
from ~ to 35% (the table in the paper by Fregeau 
et al. 2003). This is attributable to the effective dy- 
namical disruption of wide pairs and to the mergers 



of close systems. Note also that the above authors 
assumed up to 30% of the initially hard binaries in 
the GC. We choose an upper limit of 1000i? Q for 
the initial binary semiaxes, which corresponds to the 
criteria by Fregeau et al. (2003). Thus, the high (at 
first glance) percentage of the initially binary systems 
that we adopted is not overestimated and is consistent 
with available observations. Unfortunately, we can- 
not accurately estimate the number of newly forming 
binaries due to tidal captures. In any case, their per- 
centage is much lower than the fraction of the initially 
binary systems. 

The Evolution of a Binary in a Globular Cluster 

The Scenario Machine population synthesis code 
(Lipunov et al. 1996) was used to model the evolution 
of binary and single stars. We used the initial distribu- 
tions 

/(lgo) = const, (6) 

j 10.Ro ] , 
max I u \ < a < 10 3 i? Q , 

{RUM,)} 

for the semimajor axes of the binaries and 

/(Mi) oc Aff 2 - 35 , 0.1M© < Mi < 12OM (7) 

forthe mass of the primary component. Here, R^(Mi) 
is the Roche lobe radius for the primary component. 
For the component mass ratio, we used a power-law 
distribution: 

/(?)«?"', q = M 2 /M 1 <l. (8) 

Here, a q is the parameter of the distribution in com- 
ponent mass ratio. In our calculations, we varied it 
between and 1 . 

Another important parameter of the evolution of 
binaries that significantly affects the results of all 
such calculations is the energy transfer efficiency 
at the common-envelope stage, acE- Having been 
determined both by Lipunov et al. (1996) and in 
other studies of this group, it was assumed to be 
a CE — 0.5. This value corresponds to the (standard) 
value of acE = 1 in the determination used by Hurley 
et al. (2002), in which the common envelope was 
ejected only through the transfer of the gravitational 
energy of orbital motion of the approaching compo- 
nents to it. 

When describing the dynamical interactions be- 
tween cluster stars, we took into account the follow- 
ing processes in the code: 

( 1 ) The passages that led to changes in the orbital 
parameters of binaries (semimajor axis and eccentric- 

ity). 

(2) The passages that led to star exchanges. 
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(3) The passages that led to the disruption of bi- 
naries. 

(4) The interactions between single stars that led 
to the formation of binaries. 

The corresponding cross sections for these 
processes were numerically calculated by several au- 
thors (see Heggei et al. 1996; Mikkola 1984; Kim and 
Lee 1999). 

The rate of interaction between an individual bi- 
nary b at distance r from the GC center with an 
a subsystem of stars with a space density n a (r) at 
time t is given by the formula 

K ab (t) = n a (r,t)v m (r,t)a ab . 

The determination of the binary position in the GC 
and 1Z is described in more detail in the Appendix. 

To draw a change in the binding energy AE(, and 
orbital eccentricity e of a binary during the close 
passage of a third star, we used the following formulas 
for the differential interaction cross section (Heggei et 
al. 1993; Davis et al. 1992): 
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where X — AE^/E^, and the parameter A depends 
on the relative velocity of the approaching stars. The 
mean value is (X) — 0.4. 

The change in binding energy during a tripple 
collision results in momentum transfer to the binary 
barycenter. Thus, it determines the possibility of the 
binary being expelled from the cluster through such a 
collision and the distance to which the escaped binary 
can go before its merger. 

RESULTS OF THE CALCULATIONS 

Since the expulsion of merging double WDs from 
a cluster is directly related to the dynamical inter- 
actions between stars, one might expect the number 

of such systems N?^n in the evolution time to be 

J WD 2 

determined by T: 



r oc 



/ 



P0 r 'c 



P0 r ci 



where the integration is over the entire GC volume. 

The calculated numbers of merging WD pairs for 
various GC models are given in the table. We see from 
this table that the choice of a GC model (parameter T) 
most strongly affects precisely the number of WDs 

dynamically expelled from the cluster: N^ d2 oc r 65 . 
A similar power law also follows from observations for 




t, 10~ 9 yr 

Fig. 1. Merger rate of double WDs in a GC (model C, 
a q = 1) versus time, in units of 10 -6 yr -1 . The increase 
in merger rate at t ~ 7 x 10 9 yr corresponds to an in- 
crease in T during the GC core collapse: / — for systems 
merging within the GC; 2 — for systems expelled from the 
GC; and 3 — the total merger rate of double WDs in the 
GC. 



X-ray sources in the GC. In this case, the power- 
law index for X-ray sources, which are associated 
with higher-mass binaries, is close to unity: -/Vx_ ray oc 
r o.74±o.36 (p 00 iey e t al. 2003). 

The merger rate of double WDs (SN la candi- 
dates) is plotted against time for one model (C) in 



The number of merging double white dwarfs in a GC of 10 6 
stars with a total mass of 4 x 1O 5 M in an evolution time 
of 1.5 x 10 10 yr 
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381 


261 


120 
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Note. The initial number of binaries is 3 x 10 5 . The GC pa- 
rameters are given at the time of the core collapse (7 x 10 9 yr). 
The value of T was normalized in such a way that r = 1 for 
model A. iV WD 2 is the total number of merging double WDs 
in the entire evolution time; N^ 2 and N^ D2 are the numbers 
of merging double WDs within the GC and expelled from the 
GC, respectively; and v is the fraction of the expelled merging 
double WDs. For each model, the results of our calculations are 



presented for two values of a q : 
(lower rows) (see formula (8)). 



: (upper rows) and a q 
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Fig. 2. Histogram (normalized to unity) illustrating the 
distribution of WD pairs expelled from the cluster in their 
barycenter velocities (model C). 




Fig. 3. Histogram (normalized to unity) illustrating the 
distribution of WD pairs expelled from the cluster in dis- 
tance from the center of an isolated globular cluster at the 
time of their merger (model C ). 



Fig. 1. The increase in merger rate at time f~7x 
10 9 yr corresponds to a sharp increase in T during the 
evolution of the cluster (GC core collapse). 

The ratio of the total number of mergers to the 
number of evolving binaries may serve as a character- 
istic of the efficiency of the merger rate of double WDs 
that is independent of the adopted normalization (the 
total number of stars in the cluster and the fraction 
of the initially binary stars). Taking 300 mergers per 
3 x 10 5 binaries (table) as a typical value, we obtain 
an efficiency of 10~ 3 per initial binary system. This 
value is half the value obtained by Shara and Hur- 
ley (2002) (2002) in their numerical simulations and 
can be explained, in particular, by the higher efficiency 
of the common envelope (ckce — 3) in these numerical 
simulations. 

Note also that the total number of merging WDs 
as a result of the GC evolution, -/V WD 2, changes by a 
factor of only 1.5 for different GC models and is 1 00 — 
500 per cluster; most (70%) of the mergers occur after 
the collapse of the cluster core. The mean merger rate 
of double WDs in the past 5 Gyr is then ~ 10~ 13 
per year per average cluster star, which is a factor 
of ~ 3 higher than the estimated merger rate of double 
WDs in our Galaxy (Nelemans et al. 2001; Hurley et 
al. 2002). 

The distributions of WD pairs expelled from the 
cluster in their barycenter velocities and distances 
from the center of the host GC are shown in Figs. 2 
and 3. As we see from Fig. 2, the fraction of the stars 
with escape velocities of about 100 km s -1 is small (a 
few percent). 

DISCUSSION 

The Escape of Merging Double WDs from the 
System of Galactic Globular Clusters 



Our Galaxy. The GC system of our Galaxy 
extends to distances as large as 100 kpc above 
the Galactic plane and has a velocity dispersion 
of about 150 km s -1 (see, e.g., the catalogs by 
Kukarkin (1974) and Harris (1996)). The escape 
velocity from the Galactic potential is ~ 550 km s _1 
in the Galactic plane and decreases to ~ 100 km s _1 
at distances of ~ 100 kpc from the Galactic center 
(accurate estimates are difficult to obtain, because 
the form of the Galactic potential is not known at 
such distances). The escape of the stars expelled from 
the GC is preferred at the apocenter of the cluster 
orbit, near which the cluster is located for the longest 
time. Clearly, not all of the binaries expelled from the 
cluster will be able to escape from the host galaxy, but 
will remain in its halo. Thus, the total number of WD 
pairs escaping from the cluster gives an upper limit 
on the number of intergalactic type-la supernovae. 
The rough upper limit on the probable number of 
intergalactic SN la, A^sNia, that emerge during the 
mergers of double WDs expelled from the GC in a 
cluster of Af ga i galaxies is 

^SN i a = A^ WD 2 X N GC X N gal 

~ 100 x 500 x 10 3 = 5 x 10 7 

(Aqc is the number of GCs in an average galaxy). In 
a characteristic time of 10 10 yr, the formation rate of 
SNe Ia through this channel throughout the cluster 
is ~ 5 x 10~ 3 per year, which is comparable to the 
merger rate of double white dwarfs with a total mass 
exceeding the Chandrasekhar limit in a single galaxy 
(Nelemans et al. 2001). Even this upper limit is an 
order of magnitude smaller than the high fraction of 
the extragalactic SNe Ia deduced from observations. 

Dwarf elliptical galaxies. We may also assume 
that dwarf elliptical dE galaxies (the fraction of which 



THE ORIGIN OF INTERGALACTIC SUPERNOVAE 



145 



in a cluster can be significant) with a mass up to M ~ 
1O 8 M are the main source of observed intergalactic 
supernovae. The central number densities of the stars 
in such galaxies are comparable to those in dense GC 
cores. The escape velocity from a cluster at a constant 
star space density is v e oc M 1 / 3 . Therefore, applying 
our calculations to a system of 10 8 stars, we find that 
the number of escaping WD pairs from this system 
is 10 8 /10 6 x /(> v e ) x 100 - 100 (in this estimate, 
the fraction of the binaries with barycenter velocities 
higher than the escape velocity from the cluster, v e ~ 
100 km s -1 , was taken to be /(> v e ) ~ 1/100), which 
is also insufficient to explain the truly intergalactic 
supernovae. 

cD galaxies. Let us now turn to the central 
cD galaxy as the source of intergalactic supernovae. 
Although SN 1998fc in Abell 403 and SN 200 lal in 
Abell 2122/4 are projected onto the halos of the cen- 
tral cD galaxies (the projected distance is 160 kpc), 
we rule out the possibility that they belong to these 
galaxies because of the radial-velocity difference 
(Gal-Yam et al. 2003). The giant elliptical galaxies 
in the cluster centers are known to have a huge 
system of globular clusters (up to 10000). Using the 
results obtained in our calculations (~ 300 mergers 
of double WDs in cluster cores and ~ 100 merging 
pairs escaped after the core collapse), we obtain an 
average estimate of the maximum SN la rate due to 
the WD evolution in the GC in projection onto the 
halo of the cD galaxy, 400 x 10 4 /(5 x 10 9 ) - 10" 3 
per year. Only a few percent of such double WDs can 
have barycenter velocities high enough to escape from 
the galaxy; they cannot explain the observations. 

The Mergers of Double WDs in Virialized Star 
Clusters at the Centers of Galaxy Clusters 

Is there an alternative explanation for the obser- 
vations by Gal-Yam et al. (2003) that does not ap- 
peal to the evolution of 10 — 20% of the stars out- 
side galaxies? In our view, the most plausible expla- 
nation for intergalactic supernovae is the merger of 
WD pairs in a virialized system of star clusters in 
the central regions of galaxy clusters. It follows from 
observations of the central parts of the Virgo cluster 
and from simulations (see the recent review article 
by Lee (2002) and references therein) that a huge 
system of blue GCs with a large velocity dispersion 
must be formed during the formation of giant elliptical 
galaxies. Recent observations of intergalactic (up to 
250 kpc from the center) GCs in the nearby Hydra I 
and Centaurus clusters (Hilker 2002) lead to a similar 
conclusion. An enhanced density of star clusters in 
the central regions of galaxy clusters naturally arises 
during the merger of the galaxies that form giant 



cd galaxies as well as through the tidal capture of the 
GC system of the galaxies that pass near the cluster 
center. Assuming that 0.3% of the luminosity (and the 
baryon mass) in the cluster centers is concentrated in 
dense collapsed stellar structures and taking into ac- 
count the threefold increase in the merger rate of dou- 
ble WDs in such clusters, we naturally obtain ~ 1% 
of the observed SN la rate outside the visible host 
galaxies. Thus, this model predicts the existence of 
faint ((M R ) ~ -10 m ) host clusters for SN la. It is 
hoped that the purposeful search for type-la super- 
novae being conducted at present as part of various 
projects (see, e.g., Tonry et al. 2003) will improve the 
existing statistics and will reveal the faint host galax- 
ies and astrophysical sources of intergalactic SN la. 
The absence of host clusters with M R > -8 m . . . — 
10 m may serve as a critical argument for choosing 
the model of intergalactic type-la supernovae from 
intergalactic stars. 



CONCLUSIONS 

We have considered the possible formation chan- 
nel of intergalactic SN la during the mergers of dou- 
ble white dwarfs that evolve and are expelled from 
dense star clusters (GCs, the nuclei of dwarf ellip- 
tical galaxies, etc.) when they dynamically interact 
with a third star (tripple collisions). The population 
synthesis method was used to model the formation 
and evolution of close double white dwarfs in cluster 
cores by taking into account the time evolution of 
their physical parameters and to calculate the rates 
of escape of such objects from clusters. We obtained 
their distributions in escape velocities from clusters, 
in merger rates, and in distances from the host clus- 
ter after escape. We showed that the mean merger 
rate of double WDs after the core collapse could 
reach ~ 10~ 13 per year per average GC star. In a 
typical GC, the number of expelled WD pairs after 
the GC core collapse is ~ 100. For a galaxy cluster 
with several thousand members, this estimate yields 
a formation rate of intergalactic thermonuclear su- 
pernovae < 0.005 per year, which is much lower than 
the observed high rate of intergalactic supernovae of 
this type. We hypothesize that the possible observed 
high percentage of intergalactic SN la without any 
obvious host galaxies can be explained in part by 
the evolution of double WDs in dense star clusters. 
A system of such clusters naturally arises near the 
center of galaxy clusters during the formation of giant 
elliptical galaxies and, subsequently, during the tidal 
interaction between the galaxies that pass near the 
cluster center. This hypothesis can be verified through 
observations of the faint host clusters of SN la in 
nearby galaxy clusters. 
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APPENDIX 

The position of a binary in a globular cluster. 

The relative motions of stars in a GC produce 
continuous fluctuations of the gravitational field. In 
turn, these fluctuations cause the magnitude and di- 
rection of the velocity of each star to change. As a 
result, the energy E = v 2 /2 + <f>{r) and the angular 
momentum J — r x vt of the star will change. When 
calculating each evolutionary track, the changes in 
these parameters in a time interval At due to the 
gravitational interaction with the surrounding stars 
can be calculated by using the formula 



AE = 
AJ = 
At 



nei 



nji+n 1/2 y 2 ji, 
= nP(E,J), 



where y\ and y 2 , the random numbers drawn from 
a normal distribution with a mean (j/i) = (y 2 ) = 
and a variance (y 2 ) — {yi) — 1, are chosen from the 
correlation condition in such a way that 

{vm) = i 2 l^h- 

We used a method with averaging over the orbit sim- 
ilar to the method by Shapiro and Marchant (1978) 
and Marchant and Shapiro (1979) to determine the 
coefficients ei and £2,1: 
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31 



■J r„ 



(AJ}dr/v r 



3% 



(AJ 2 )dr/v r 



= 2 f 

-- 2 f 3 {AEA J)dr/v r . 



coefficients: 
(AE) 



viAv^ + ^AvD + ^Av 2 ), 



(AE 2 )=v 2 (Av 2 ), 



J 



(AJ) = —(A 
J 2 



AJ 



^ J ' Z ) = + 2 {r 2 -^)(A«i), 

(AEAJ) = J(Avf { ) 

For a Maxwellian velocity distribution of the field 
stars, the diffusion coefficients reduce to standard 
form(Spitzer 1987): 



1 + 



rn 



(A V]] ) = -2 

((A V]] ) 2 
((Av ± ) 2 ) = 2n f Tk 



rrifj 
2n f Tk 



) n f Tk 2 G(x), 
G{x) 



$(x) - G{x) 



where <&(x) is the error function 

X 

2 



and 



G(x) 



r l/2 



/ 



exp 



-y 2 dy 



<&(x) - x$'(x) 
2x~ 2 ' 



Using E and J, we can determine the orbital pa- 
rameters of the star within the GC: the orbital period 



P(E,J) = 2 



/ dr/v r , 

J rp 



the pericenter r p and the apocenter r a of the stellar 
orbit in the cluster from the condition 

v 2 = 2E- 24>{r) - J 2 /r 2 = 0. 



The rate of interaction between a binary b and an 
a subsystem of stars at time t is given by the formula 

fc a b{t) = x / n a (r,t)v(r,t)/v r dr. 

r{&,J) Jrp 



Here, r p and r a are the pericenter and apocenter of 
the stellar orbit in the cluster. These coefficients are 
the mean rates of change in E and J due to several 
approaches over the orbital period and the mean cu- 
mulative values over the orbital period. 



REFERENCES 

1. M. Bellazzinni, P. F. Fusi, P. Montegriffo, et al., As- 
tron. J. 123,2541 (2002). 

2. E. Capellaro, M. Turatto, D. Yu.Tsvetkov, et al., As- 
tron. Astrophys. 322, 431 (1997). 



THE ORIGIN OF INTERGALACTIC SUPERNOVAE 



147 



3. N. N. Chugai and L. R. Yungel'son, Pis'ma Astron. 
Zh. (2004, in press); astro-ph/0308297 (2003). 

4. N. V. Dunina-Barkovskaya, V. S. Imshennik, 
S. I. Blinnikov, Pis'ma Astron. Zh. 27, 412 (2001) 
[Astron. Lett. 27, 353(2001)]. 

5. M. B. Davis, W. Benz, and J. G. Hills, Astrophys. J. 
424,870(1992). 

6. J. M. Fregeau, M. A. Gurkan, K. J. Joshi, 
and F. A. Rasio, Astrophys. J. (in press); astro- 
ph/0301521 (2003). 

7. M. Hamyu, M. M. Phillips, N. B. Suntzeff, et al, 
Nature (in press); astro-ph/0306270 (2003). 

8. W. E. Harris, Astrophys. J. 112, 1487 (1996). 

9. D. C. Heggei, C. Douglas, and P. Hut, Astrophys. 
J. Subbpl. Sen 85, 347(1993). 

10. D. C. Heggie, P. Hut, and S. L.W. McMillan, Astro- 
phys. J. 467, 359 ( 1996). 

1 1 . M. Hilker, astro-ph/02 1 0566 (2002). 

12. F. Hoyle and W. A. Fauler, Astrophys. J. 132, 565 
(1960). 

13. J. R. Hurley, C. A. Tout, and O. R. Pols, Mon. Not. R. 
Astron. Soc. 329, 897(2002). 

14. I. Iben, Jr. and A. V. Tutukov, Astrophys. J., Suppl. 
Ser. 54,355(1984). 

15. A. Gal-Yam, D. Maoz, P. Guhathakurta, and 
A. V. Filippenko, Astron. J. 125, 1087 (2003). 

16. S.S.Kim and H.M.Lee, Astron. Astrophys. 347, 123 
(1999). 

17. S. S. Kim, H. M. Lee, and J. Goodman, Astrophys. J. 
495,768(1998). 

18. B. V. Kukarkin, Globular Star Clusters (Nauka, 
Moscow, 1974) [in Russian]. 

19. M. G. Lee, J. Korean Astron. Soc. 35, 1 (2002). 

20. W. Lee, A. V. Filippenko, R. R. Treffers, et al., Astro- 
phys. J. 546, 734(2001). 

21. V. M. Lipunov, K A. Postnov, and M. E. Prokhorov, 
Astrophys. Space Sci. Rev. 9, 1 (1996). 

22. M. Livio and A. Riess, submitted to the Astrophys. J.; 
astro-ph/0308018. 

23. A. B. Marchant and S. L. Shapiro, Astrophys. J. 234, 
317(1979). 



24. R. W. Michie, Mon. Not. R. Astron. Soc. 125, 127 
(1063). 

25. S. Mikkola, Mon. Not. R. Astron. Soc. 208, 75 
(1984). 

26. R. Napiwotzki, N. Christlieb, H. Drechsel, et al., 
Messenger 112,25(2003). 

27. G. Nelemans, L. R. Yungelson, S. F. Portegies Zwart, 
and F. Verbunt, Astron. Astrophys. 365, 491 (2001). 

28. K. Nomoto, Astrophys. J. 253, 798 (1982). 

29. S. Perlmutter, G. Albering, G. Goldhaber, et al., As- 
trophys. J. 517, 565(1999). 

30. D. Pooley, W. H. Lewin, S. F. Anderson, et al., Astro- 
phys. J. 591, 131 (2003). 

31. A. G. Riess, A. V. Filippenko, P. Challis, et al., As- 
tron. J. 116, 1009(1998). 

32. E. H. Rubinstein and C. D. Bailyn, Astrophys. J. 474, 
701 (1997). 

33. S. L. Shapiro and A. B. Marchant, Astrophys. J. 225, 
603(1978). 

34. M. M. Shara and J. R. Hurley, Astrophys. J. 571, 830 

(2002) . 

35. E. I. Sorokina, S. I. Blinnikov, O. S. Bartunov, Pis'ma 
Astron. Zh. 26, 90 (2000) [Astron. Lett. 26, 67 
(2000)]. 

36. L. Spitzer, Jr., Dynamical Evolution of Globular 
Clusters (Princeton Univ. Press, Princeton, 1987; 
Mir, Moscow, 1990). 

37. J. L. Tonry, B. P. Schmidt, B. Barris, et al., Astrophys. 
J. (in press); astro-ph/0305008 (2003). 

38. N. Trentham and R. B. Tully, Mon. Not. R. Astron. 
Soc. 335, 712(2002). 

39. S. van den Berg, Li Weidong, and A. Fillippenko, 
Publ. Astr. Soc. Pac. (in press); astro-ph/0308195 

(2003) . 

40. R. F. Webbink, Astrophys. J. 227, 355 ( 1984). 

41. J. Whelan and I. Iben, Jr., Astrophys. J. 186, 1007 
(1073). 



4 



